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1. Introduction 



The finite temperature transition for QCD at zero baryon density has been a subject of intense 
study for the lattice gauge theory community. The nature of the phase transition and the value for 
the pseudo-critical temperature T c is not only of fundamental interest, but also has phenomenologi- 
cal implications for experimental studies of relativistic heavy ion collisions, as well as the evolution 
of the early universe. 

Most large-scale studies of the QCD transition on the lattice have used staggered fermions. 
Staggered fermions are computationally inexpensive, and preserve a U(l) remnant of the SU (2)l<S> 
SU (2)r chiral symmetry at finite lattice spacing. However, lattice artifacts in the staggered formu- 
lation break flavor symmetry, leading to, e.g. artificial mass splittings in the hadron spectrum. 

It would be desirable to use a lattice action that more faithfully reproduces the continuum 
symmetries of QCD. One such action is known as domain wall fermions (DWF). DWF possesses 
an exact SU (2) g) SU (2) chiral symmetry, even at finite lattice spacing. This is accomplished by 
adding a fictitious fifth dimension into the lattice action, where the physical chiral modes are bound 
to opposite four-dimensional walls. If the fifth dimension is infinitely large, the chiral modes are 
isolated and the chiral symmetry is exact. However, for any real lattice calculation, one must choose 
a finite fifth-dimensional size, L s . In such a calculation, L s controls the amount of mixing between 
the chiral modes, and thus the residual chiral symmetry breaking. 

There have been some past studies of chiral symmetry breaking with domain wall fermions 
|2]]. However, in these past studies, it has been difficult to control the residual chiral symmetry 
breaking, which increases very quickly as one moves towards strong coupling. In this work, we 
present a new study of the transition region with domain wall fermions, employing a modified 
gauge action, the Dislocation Suppressing Determinant Ratio (DSDR) with the Iwasaki action on 
the gauge links. This combination keeps residual chiral symmetry breaking small, allowing us to 
study the transition region while keeping the physical pion mass fixed, m % rj 200 MeV. 



2. Method 

2.1 Dislocation Suppressing Determinant Ratio (DSDR) 

In our lattice calculation, we use domain wall fermions, along with the DSDR method and 
the Iwasaki gauge action. The finite fifth dimension results in a residual chiral symmetry breaking 
caused by a mixing of the chiral modes between the four-dimensional walls. To leading order, 
this results in a simple additive renormalization to the quark masses, the residual mass wires- The 
residual mass may be parameterized by Q : 

e —AcL s J 

^res =c\p H {X c )— \-c 2 p H {Qi) — , (2.1) 

As As 

where p#(A) denotes the eigenvalue density of the effective four-dimensional Hamiltonian. The 



first term in eqn. 2.1 comes from de-localized states with eigenvalues near the mobility edge, X c , 
while the second term comes from localized "lattice dislocations" with near zero eigenvalue. As 
one moves to strong coupling, the rapid proliferation of these localized dislocations means that 



the latter term in eqn. 2.1 dominates the residual mass for moderate values of L s ~ 16 — 32. In 
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this regime, m r es ~ l/L s , meaning that it becomes very inefficient to reduce the chiral symmetry 
breaking by increasing the size of the fifth dimension. 

The near-zero eigenmodes of the effective Hamiltonian are closely related to the near-zero 
eigenmodes of the the Hermitian Wilson-Dirac operator, H$ = y 5 Dvi/(— Mo), where Mo is the do- 
main wall height. Thus, it has been suggested that augmenting the action with the determinant of 
H 5 will suppress exactly those modes that make the biggest contribution to m res [Q, ||]. However, it 
is exactly these near-zero modes of H5 that allow topology to change during a molecular dynamics 
evolution. In order to allow for topology change during our Monte Carlo calculation, we instead 
introduce the following weight factor into our action: 



W(M ,s h ,Sf) 



det 



D^i-Mo + ie f f)D w (-Mo + iEff) 



det 



D^i-Mo + i£ b r 5 )D w (-M + ie b y 5 ) 



(2.2) 



det 


Dl(-M )D W ( 


-M ) + ej 


det 


dI(-m )d w ( 


-Mo) + e 2 b 



where e b and Ef are chirally-twisted mass terms. By tuning e b and Ef appropriately, we can find 
a weighting factor that reduces ra r es sufficiently without completely eliminating topology change 
during our HMC evolution. The choice of Ef = 0.02 and £ b = 0.50 satisfies these conditions [||]. We 
call the weighting factor, W(Mo,£b,£f) the Dislocation Suppressing Determinant Ratio (DSDR). 

2.2 Line of constant physics 

Using the DSDR method, we can choose our input bare light quark masses mi such that the 
total quark mass m[ [ = mi + m re s corresponds to a fixed physical pion mass, m n k, 200 MeV. We 
have produced configurations at 7 different temperatures with spatial size N a = 16 and temporal 
extent N x = 8, spanning the temperature range T G [139, 195] MeV. This corresponds to lattice 
spacings of a « 0. 13 — 0. 18 fm. For each ensemble, we have produced between 2996-7000 molec- 
ular dynamics trajectories. Table |] summarizes the parameters that we use in our finite temperature 
ensembles. 

Fig. [I] shows a comparison of m r es as a function of temperature between our current calcula- 
tion with domain wall fermions and the DSDR method, compared to domain wall fermions without 
the DSDR method, both at N r = 8. With L s = 32, the DSDR method offers a substantial improve- 
ment in residual chiral symmetry breaking. In fact, for temperatures T > 170 MeV, one needs to 
have at least L s = 96 in order to achieve residual masses equal to what can be obtained using DSDR 
at L s = 32. 

In order to determine the lattice scales in our calculation and the quark masses needed in our 
finite temperature ensembles, we have also used results from zero temperature calculations at three 
values of the gauge coupling /3 = 1.70,1.75,1.82 that lie within our range of interest. Table || 
summarizes the parameters for our zero temperature ensembles. 

Using interpolations and extrapolations from the zero temperature data, we are able to deter- 
mine the total bare light quark masses, rri[ [ = mi + m re s required to keep the physical pion mass 
fixed to m n = 200 MeV. Fig. |2| shows this value of m[ Q [ as a function of the bare coupling, /3. 
Also shown are the values of m re s obtained on the DSDR finite temperature lattices with L s = 32 



3 



The finite temperature phase transition from domain wall fermions 



Michael Cheng 



and L s = 48. Even with L s = 48, the lowest temperature ensemble (/3 = 1.633) requires a negative 
input light quark mass. At j8 = 1.671 we have generated two separate ensembles with L s = 32 and 
L s = 48, with a negative and positive input light quark mass, respectively, in order to test if the neg- 
ative quark mass produces any deleterious effects. With these parameters, we found no evidence 
of any "exceptional configurations" that are, in principal, possible with a negative quark mass. 



Figure 1: Comparison of mres from domain 
wall fermions without the DSDR method at L s = 
32,96 and what is obtained using DSDR atL v = 
32. 




Figure 2: The solid black curve plots the total 
bare quark mass m[ [ required for a fixed physi- 
cal pion mass m K = 200 MeV as a function of the 
bare coupling j3 . The red curve and points show 
'«res for the L s = 32 DSDR lattices. The blue 
points show mres for the L s = 48 ensembles. 



T (MeV) 


P 


Ls 


m\a 


Finite Temperature Ensembles 
m s a m r es m-n (MeV) 


Xl,disc/T 2 


V MS l T 2 
M. disci 1 


Traj. 


139(6) 


1.633 


48 


-0.00136 


0.0519 


0.00588(39) 


191(7) 


37(3) 


17.2(14) 


2996 


149(5) 


1.671 


32 


-0.00189 


0.0464 


0.00643(9) 


199(5) 


44(3) 


19.9(10) 


6000 


149(5) 


1.671 


48 


0.00173 


0.0500 


0.00295(3) 


202(5) 


41(2) 


18.5(9) 


7000 


159(4) 


1.707 


32 


0.000551 


0.0449 


0.00377(11) 


202(3) 


43(4) 


18.8(18) 


3659 


168(4) 


1.740 


32 


0.00175 


0.0427 


0.00209(9) 


197(2) 


35(5) 


14.9(21) 


3343 


177(4) 


1.771 


32 


0.00232 


0.0403 


0.00132(6) 


198(2) 


25(4) 


10.4(17) 


3540 


186(5) 


1.801 


32 


0.00258 


0.0379 


0.00076(3) 


195(3) 


11(4) 


4.5(1.6) 


4715 


195(6) 


1.829 


32 


0.00265 


0.0357 


0.00047(1) 


194(4) 


5(3) 


2.0(1.2) 


6991 



Table 1: Summary of finite temperature ensembles with DSDR and results for the disconnected chiral 
susceptibility, Xl,disc- 



3. Chiral Symmetry Restoration 

On our finite temperature lattice ensembles we have computed the connected, spatial two- 
point functions in various mesonic channels, Gr(x). From these, we can trivially calculate the 
corresponding connected susceptibilities, 

ycon 

^ r =iV T 2 £Gr(x). (3.1) 
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Zero Temperature Ensembles 



p 


No 


N? 


L s 


mia 


m s a 


"Jres 


m n (MeV) 


ro/a 


a~ l (GeV) 


Traj. 


1.70 


16 


32 


32 


0.013 


0.047 


0.00420(2) 


394(9) 


2.895(11) 


1.27(4) 


1360 


1.70 


16 


32 


32 


0.006 


0.047 


0.00408(6) 


303(7) 


2.992(27) 


1.27(4) 


1200 


1.75* 


32 


64 


32 


0.0042 


0.045 


0.00180(5) 


246(5) 


3.349(20) 


1.36(3) 


1288 


1.75* 


32 


64 


32 


0.001 


0.045 


0.00180(5) 


172(4) 


3.356(22) 


1.36(3) 


1560 


1.82 


16 


32 


32 


0.013 


0.040 


0.00062(2) 


398(9) 


3.743(28) 


1.55(5) 


2235 


1.82 


16 


32 


32 


0.007 


0.040 


0.00063(2) 


304(7) 


3.779(37) 


1.55(5) 


2134 



Table 2: Summary of zero temperature ensembles with DSDR. Lattice scales determined using ro = 
0.487(9) fm, after extrapolation to the chiral limit. *: /J = 1.75 results are RBC-UKQCD preliminary 
results. 

We have also measured the one-flavor scalar chiral condensate, (Yi¥i) an d the pseudoscalar 
condensate, (yf// V/) with a stochastic estimator. Using these, we can compute the disconnected 
scalar and pseudoscalar susceptibilities, 



v disc 

M 

T 2 

.disc 



A^X 3 (<(¥/W) 2 >-(¥/W) 2 ) (3-2) 
Xr -NlNU((W^Yi) 2 )-(W,f¥i) 2 ) (3-3) 



T 2 

With the connected and disconnected susceptibilities, we can trivially reconstruct the flavor 
singlet (a) and non-singlet (5) susceptibilities in the scalar channel, 

nj v con I v disc v v con 

j2~ j2 ' j2 ~ J2 ' ' 

and also the flavor singlet (T]) and non-singlet {%) susceptibilities in the pseudoscalar channel, 

v v con v disc v v con 

AT) _ A5 A5 Aft _ /T 

J2~ j2 ' j2 ~ j2 ' > 

One of the signatures for chiral symmetry restoration is the divergence of the disconnected 
chiral susceptibility (xf' sc ) at T c , along with the corresponding 0(4) scaling near T c , in the chiral 
limit. Thus, one method to identify the crossover temperature at finite quark mass is by the peak 
in the chiral susceptibility. Fig. [3] is a comparison of the disconnected chiral susceptibility in our 
current calculation with the earlier results from which were not performed along a line of 
constant physics. Although there is a clear peak near T m 170 MeV in the earlier results, the results 
in the low temperature region are significantly distorted compared to our current calculation. In our 
calculation, a peak is visible near T « 160 MeV. 

Fig. [| compares our results with the results of various improved staggered actions 0. In 
order to compare results from different actions, all results are converted into the MS renormaliza- 
tion scheme. While the lightest pions in the staggered calculation correspond to a physical quark 
mass that is approximately half as light as in our calculation, there seems to be reasonably good 
agreement for Xi f° r T > T c . However, whereas the location of the peak in %f KC , T c w 160 MeV 
is approximately the same as for the HISQ action with N r = 12, the magnitude of the susceptibility 
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V DWFDSDR 

□ Asqtad N =8 

■ Asqtad N =12 

* HISQ N =6 
O HISQN = 

• HISQ N =12 




180 

T (MeV) 



Figure 3: Comparison of the disconnected chiral 
susceptibility in our calculation (DWF DSDR) 
with earlier results (DWF L s = 32), where 
the finite temperature ensembles do not lie along 
a line of constant physics. 



170 180 
T (MeV) 



210 220 



Figure 4: Comparison of the disconnected chi- 
ral susceptibility with results from various stag- 
gered actions [Q]. All results have been con- 
verted to the MS renormalization scheme. 



for T < T c is significantly greater for DWF than the various staggered actions. This may be due 
to finite volume effects, as the DWF calculation has an an aspect raio Na/N^ = 2, compared to 
N a /N x = 4 for the staggered calculations. 

In the chirally -restored phase, we can also derive relations between the scalar and pseudoscalar 
susceptibilities. The chiral transformation mix the scalar and pseudoscalar channels, so that in the 
high temperature phase, we have the relations % n = % a and also Xs = Xn • Thus, we expect that the 
differences in these susceptibilities should vanish in the chiral limit at temperatures where chiral 
symmetry is restored. 



Using the definitions in eqns. 3.4 and fO| , 
we can also show that the difference in the 
pseudoscalar and scalar disconnected sus- 
ceptibilities, Adisc = %f sc — %f lsc should also 
vanish in the chirally restored phase. Fig. 
|5| shows Adis C /T 2 as well as the difference 
(Xn —%<j)/T 2 - Even though we are not in 
the chiral limit, both of these differences are 
already consistent with zero for T pa 170 — 
1 80 MeV, suggesting that chiral symmetry is 
already well restored at those temperatures. 



180 
160 
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100 
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60 
40 
20 




T[MeV1 ^_ 

140 150 160 



Figure 5: Plot of A c i isc /T 2 = {%f sc - xf sc ) /T 2 
and (Xk — Xo)/T 2 as a function of temperature. 
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4. Conclusion 

In this work, we have examined the restoration of chiral symmetry in finite temperature QCD 
using domain wall fermions on lattices of spatial size N a = 16 and temporal extent 7V T = 8 in the 
temperature range T G [139, 195] MeV. In order to control the residual chiral symmetry breaking, 
we have employed the DSDR method, augmenting the gauge action with a weighting factor that 
suppresses the localized dislocations that contribute most to m r es- By controlling residual chiral 
symmetry breaking, we were able to use a fixed, physical pion mass m n pa 200 MeV. 

In order to find the chirally restored phase, we examined the light disconnected chiral sus- 
ceptibility, which exhibits a peak at the crossover temperature T c . A comparison with an earlier 
domain wall calculation, which did not have a fixed pion mass, show the distortions introduced if 
one does not do calculations along a line of constant physics. A comparison with recent results with 
improved staggered fermions shows significant differences in the chiral susceptibility, especially 
for T < T c . Another signature of chiral symmetry restoration that we examined were the scalar 
and pseudoscalar susceptibilities, both the flavor singlet and non-singlet channels. By examining 
the differences in these susceptibilities, we could see apparent chiral restoration for T pa 170 — 180 
MeV. 

In addition to the restoration of chiral symmetry, we have also examined the effective restora- 
tion of U{\)a symmetry as well as the eigenvalue spectrum on these finite temperature lattices. 
These results are presented in these proceedings in [fy respectively. Future calculations are 
underway to extend these results to larger spatial volume in order to control finite-volume effects, 
as well as to smaller quark masses, so that we can examine the quark mass dependence of these 
quantities. 

References 

[1] P. Chen, et. al., The Finite temperature QCD phase transition with domain wall fermions, Phys. Rev. 
D64014503 [hep-lat/0006010]. 

[2] M. Cheng, et. al., The finite temperature QCD using 2 + 1 flavors of domain wall fermions at N z = 8, 
Phys. Rev. D81 054510 [arXiv: 0911 . 3450]. 

[3] D.J. Antonio, et. al, Localization and chiral symmetry in three flavor domain wall QCD, Phys. Rev. 
D77 014509 [arXiv: 0705.2340]. 

[4] P. Vranas, Gap Domain Wall Fermions, Phys. Rev. D74 034512 [hep-lat / 6 6 1 4]. 

[5] H. Fukaya, et. al., Lattice gauge action suppressing near-zero modes ofH(W), Phys. Rev. D74 094505 

[hep-lat/060702]. 

[6] D. Renfrew, Controlling Residual Chiral Symmetry Breaking in Domain Wall Fermion Simulations, in 
proceedings 

[7] A. Bazavov, et. al., The chiral and deconfinement aspects of the QCD transition, 



[arXiv: 1111.1710]. of Lattice 2008, PoS (LATTICE2008) 04 8 [arXiv: 0902.2587] 



[8] P. Hegde, Looking for U(1)a Restoration in Hot QCD with Domain Wall Fermions, in proceedings of 



Lattice2011, PoS (Lattice 2011) 014 [arXiv: 1112 . 0364] 



[9] Z. Lin, Dirac Eigenvalue Spectrum at Finite Temperature Using Domain Wall Fermions, in 



proceedings of Lattice 2011, PoS (Lattice 2011)204 [arXiv :1111.0988] 



7 



